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ABSTRACT 



Context, BA-type supergiants show a high potential as versatile indicators for modern astronomy. This paper constitutes the first in a 

■ series that aims at a systematic spectroscopic study of Galactic BA-type supergiants. Various problems will be addressed, including 
in particular observational constraints on the evolution of massive stars and a determination of abundance gradients in the Milky Way. 
Aims. The focus here is on the determination of accurate and precise atmospheric parameters for a sample of Galactic BA-type su- 
pergiants as prerequisite for all further analysis. Some first applications include a recalibration of functional relationships between 
spectral-type, intrinsic colours, bolometric corrections and effective temperature, and an exploration of the reddening-free Johnson Q 
and Stromgren [ci] and /3-indices as photometric indicators for effective temperatures and gravities of BA-type supergiants. 
Methods. An extensive grid of theoretical spectra is computed based on a hybrid non-LTE approach, covering the relevant parameter 
space in effective temperature, surface gravity, helium abundance, microturbulence and elemental abundances. The atmospheric pa- 
rameters are derived spectroscopically by line-profile fits of our theoretical models to high-resolution and high-S/N spectra obtained 
at various observatories. Ionization equilibria of multiple metals and the Stark-broadened hydrogen and the neutral helium lines con- 
stitute our primary indicators for the parameter determination, supplemented by (spectro-)photometry from the UV to the near-IR. 
Results. We obtain accurate atmospheric parameters for 35 sample supergiants from a homogeneous analysis. Data on effective tem- 
peratures, surface gravities, helium abundances, microturbulence, macroturbulence and rotational velocities are presented. The inter- 
stellar reddening and the ratio of total-to-selective extinction towards the stars are determined. Our empirical spectral-type-r cff scale 
is steeper than reference relations from the literature, the stars are significantly bluer than usually assumed, and bolometric corrections 
differ significantly from established literature values. Photometric r c(f -determinations based on the reddening-free (9-index are found 
to be of limited use for studies of BA-type supergiants because of large errors of typically ±5% (lcr statistical)±3% (lcr systematic), 
compared to a spectroscopically achieved precision of 1-2% (combined statistical and systematic uncertainty with our methodology). 
The reddening-free [ci ]-index and /3 on the other hand are found to provide useful starting values for high-precision/accuracy analyses, 
with uncertainties of ±1% ±2.5% in T sff , and ±0.04±0. 13 dex in logg (lcr-statistical, lcr-systematic, respectively). 

Key words. Stars: atmospheres - Stars: early-type - Stars: fundamental parameters - Stars: rotation - supergiants - dust, extinction 
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■ 1. Introduction This can provide observational constraints on galactochemical 

evolution models of different galaxy types as well as on stel- 
O Supergiants of late B and early A-type (BA-type supergiants) lar evolution models for a wide range of met allicities, see e.g. 
CN ■ are among the visually brightest stars, reaching absolute magm- Przvbilla (20 08) for a review, and in a broader context on the 
r-i ■ tudes of up to My » -9.5. Therefore, they show a high potential galaxy mass _ me tallicity relationship ( KudritzkiitiD|20i2). 
> . as versatile indicators for stellar and galactic studies oyer large Moreover> B A-type supergiants can act as standard candles 
distances using ground-based telescopes of the 8-10m class. The for difjtance determinations via use of the wind momentum- 
focus of quantitative studies of individual BA-type supergiants luminosit relationship (Puis et al. 1996; Kudritzki et al. 1999) 
H . lies in extragalactic research at present. Objects in many of the orjieflux-weighted gravity-luminosity relationship (FGLR, 
star-forming galaxies of the Local Group and even beyond have Kudritzki et al 2003 200 8b). The latter is of particular inter- 
been investigated at hi gh and intermed i ate spectral resolution ^ as ^ steUar metallicit and me inters tellar reddening are 



X 



(see e. S the rev iews by Venn et alj2003 |Kudritzki et al. || 2008at also derived fa me antitative anal is _ In conseqU ence, sources 

| Kudntzki||2010 and references therein). Such observations in of tematic errors that trouble photom etric indicators like the 

other galaxies allow the spatial distribution of elemental abun- Cepheid p eriod . lum inosity relationship (Freedman et al. 2001; 

dances to be mapped, using individual supergiants as tracers. |Kudritzki et alJl2Tjrj8bl) can be constrained by this spectroscopic 

method. Systematic studies therefore promise refinements to 

* Based on observations collected at the Centra AstronomicoHispa- ^ extragalactic distance scale and to the determinatio n of the 

no Aleman at Calar Alto (CAHA), operated lointly by the Max-Planck TT , , . , , i T , , .. , . „ TT , ni-irvnk 

y . <-.. i . j., t .... , . . n, ■ j a j i - Hubble constant to be achieved Kudritzki & Urbaneia 2012). 

Institut lur Astronomie and the Instituto de Astrohsica de Andalucia " 

(CSIC), proposals H2001-2.2-01 1 and H2005-2.2-016. Accordingly, detailed quantitative studies of nearby stars of 

** Based on observations obtained at the European Southern Obser- mis kind > in the Milk y Wa Y' should be of 8 reat interest. Modern 

vatory, proposals 62.H-0176 and 079.B-0856(A). Additional data were spectrographs open up the possibility to obtain spectra of excep- 

adopted from the UVES Paranal Observatory Project (ESO DDT tional quality of these bright stars, even with small telescopes 

Program ID 266.D-5655). of the 2 m-class. Galactic targets may therefore act as testbeds 
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for scrutinising stellar atmosphere analysis techniques for BA- 
type supergiants before using them in extragalactic applications 
dPrzv billa 2004 iPrzvbilla et al .1120061) . However, there is much 
more to gain from a systematic study of BA-type supergiants in 
the Milky Way. New light may also be shed on some facets of 
the actions of the cosmic cycle of matter in the high-metallicity 
environment of a typical giant spiral galaxy. 

Studies of stellar structure and evolution are the basis for 
our understanding of the synthesis of the heavy elements (metal 
yields) and their injection into the cosmic matter cycle by stel- 
lar winds and supernovae. The most recent generations of evo- 
lution models for massive stars that account for effects of rota- 
tion and mass loss as well as magnetic fields (see e.g. the re- 
cent review by Ma eder & Mevnetll2012l and references therein) 
are highly successful in describing many observational aspects 
of the massive star populations in general. However, many de- 
tails - in particular related to the subtle signatures of mixing of 
CNO-cycled produ cts in the stars - ar e subject to intense d ebate 
at present, see e.g.lHunter et al.1 d2009t) . lMaeder et al.1 (120091) and 
IPrzvbilla et al.l d2010l) . Here, highly accurate observational data 
on fundamental stellar parameters and light element abundances 
are required for thoroughly testing the models. BA-type super- 
giants are primary candidates for the study of the post-main- 
sequence evolution, linking core H-burning main-sequence stars 
with core He-burning red supergiants (some He-burning objects 
may also be on a blue loop). 

Furthermore, studies of stars are the key to understand the 
chemical evolution of the Galaxy. Luminous BA-type super- 
giants are highly useful in mapping elemental abundance pat- 
terns throughout the Milky Way. Systematic investigations of 
these short-lived stars can constrain the present-day end point 
of ~13Gyr of Galactic evolution in form of abundance gra- 
dients throughout the Galactic disk. These represent one of 
the major observational constraints to Galactoche mical evolu- 
tion models. Previous work using H n-regions (e . g . IShaver etafl 
119831: lEsteban et all 120051: iRudolph et all 12006. and references 
therein), B-type main sequence stars (e.g. Gummersbach et all 
[19981: iRolleston et alj|2000l:lpaflon & Cunh a 2004) or Cepheids 
(e.g. lAndrievskv et al.l2004clPedicellie t al. 2009, and references 
therein) may be verified, complemented and extended by studies 
of BA-type supergiants. Precise and unbiased elemental abun- 
dances for a larger sample of objects with well-constrained dis- 
tances, distributed throughout the disk, are required for this. 

While the literature on properties of hot, massive stars of 
O- and early B-type is rich, few studies have concentrated on 
tepid BA-type supergiants in the past ~20 years, where model 
atmosphere techniques have reached some kind of maturity. 
Pioneering work was do ne bv lVennl d!995albl) . later updated by 
IVenn & Przvbillal (|2003), who determined atmospheric param- 
eters and chemical abundances for a sample of 22 Galactic A- 
type supe rgiants, mostly l ess-lum inous ones of luminosity class 
II and lb. IVerdugo etall (i 1999b) derived basic stellar parame- 
ters for 31 early A-type supergiants, while the compilation of 
iLvubimkov et al.l yOlO) provides such data for 8 mostly less- 
luminous late A-type supergiants. Moreover, atmospheric pa- 
rameters for a smaller numbe r of lat e B-t ype supergiants wer e 
provided by iMcErlean et all d 1999b and iFraser et alj ((2010). 
Finally, the work of iTakeda & Ta kada-Hidail (120001 and refer- 
ences therein) on a representative sample of BA-type supergiants 
is the most comprehensive one concerning abundance determi- 
nations that account for deviations from the standard assumption 
of local thermodynamic equilibrium (non-LTE) so far. 

However, a comprehensive and homogeneous study of 
Galactic BA-type supergiants in which atmospheric and funda- 



mental stellar parameters, and abundances for light, a-process 
and iron group elements alike are derived is still unavailable. 
This was our motivation to embark on the present work, com- 
bining state-of-the-art observational data based on echelle spec- 
tra with a sophisticated hybrid non-LTE analysis methodol- 
ogy dPrzvbilla et al J 120061 with extensions to facilitate a semi- 
automatic analysis as described below). The aim is to provide 
as accurate (i.e. minimised systematic uncertainties) and pre- 
cise (i.e. minimised statistical errors) data as presently possible 
for a larger number of objects. Here, in Paper I, we introduce 
the observational database and derive the atmospheric parame- 
ters of the sample stars. This constitutes the basis for all further 
discussion, that will concentrate on observational constraints to 
evolution models for massive stars (Paper II) and to models of 
Galactochemical evolution (Paper III). 

The paper is organised as follows. In Sect. 2 we describe the 
available observational data, while our methods for the determi- 
nation of stellar parameters and elemental abundances are dis- 
cussed in Sect. 3. We present the results of our studies in Sect. 4 
and some first applications in Sect. 5. Finally, a short summary 
is given in Sect. 6. 

2. Observations and data reduction 

While high-resolution and high signal-to-noise (S/N) spectra are 
required in order to achieve the high accuracy and precision we 
aim for, a near-complete wavelength coverage from ~3900 to 
9100 A is also essential for our analysis. Only then all the strate- 
gic sets of lines both for the parameter and the abundance anal- 
ysis are available in all the stars. 

We obtained data on 35 bright stars (V < 8.7 mag) using three 
echelle spectrographs both on the northern and southern hemi- 
sphere, thus sampling a wide range in Galactic longitude. The 
targets were chosen in a way that covers the examined parame- 
ter domain (B8 to A3 in spectral type, lb to la in luminosity class, 
according to an initial classification from the SIMBAD database 
at CDSlT) rather homogeneously. However, overall more lumi- 
nous objects were favoured in order to facilitate larger distances 
to be reached within our study. Stars in both OB associations 
(25 objects) and in the field (10 objects) were observed. The 
resulting list of program stars, their spectral classification, as- 
sociation membership, photometric properties and observational 
details (observation dates, exposure times and S/N-ratio of the 
spectra) a re summarised in T able Q] Note that objects in com- 
mon with lPrzvbilla et alJ d2006l) and lschiller & Przvbillal {2008) 
were reanalysed for the sake of homogeneity (see that papers on 
details of the data reduction). 

Most of the spectra were obtained in observing runs 
with the Fibre Opt ics Echelle Cassegrain Spectrograph (Foces, 
|Pfeifferetai1 [i998) on the Calar Alto 2.2 m telescope in 2001 
and 2005. These cover a wavelength range from 3860 to 9580 A 
at a resolving power R-A/AA^40000. Only the data for 
HD 195324 were acquired with a different setup, trading higher 
resolution for a lower wavelength coverage. Relatively bright ob- 
jects were observed in order to reach the desired spectrum qual- 
ity of more than 150 in S/N ratio. 

A semi- automatic pipeline w as used for the reduction of the 
Foces data (Pfeiffer et al.|[T998l) . performing subtraction of bias 
and dark current, flatfielding, wavelength calibration, rectifica- 
tion and merging of echelle orders. As the Foces observations 
were obtained over several nights in different years, the condi- 
tions and quality varied for each set of exposures. Before the data 



http : // simbad . u-strasbg . fr/ simbad/ 
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Table 1. The star sample: id, spectral type, OB association membership" , photometry* , and observational details c . 



# 


Object 


Sp.T" 


Sp.T" 


OB Assoc. 


V 


B- V 


U -B 


Date 


f 


S/N v 












mag 


mag 


mag 




s 




Foces 


R = 40 000 




















1 


HD 12301 


B8I1/ 


B8Ib 


Field 


5.589±0.011 


0.370±0.014 


-0.276+0.009 


30/09/2001 


480 


203 


2 


HD 12953 


Al Iae-f 


Allae 


Per OBI 


5.691±0.021 


0.614±0.007 


-0.014±0.009 


26/09/2001 


300 


230 


3 


HD 13476 


A3 lab 


A3 lab 


Per OBI 


6.431+0.020 


0.600+0.013 


0.220±0.028 


30/09/2001 


900 


202 


4 


HD 13744 


AO lab 


AO lab 


Per OBI 


7.592±0.014 


0.741±0.012 


0.180±0.000 


27/09/2005 


2700 


182 


5 


HD14433 


AIL/ 


Alia 


Per OBI 


6.401±0.019 


0.567±0.008 


0.030±0.010 


30/09/2001 


600 


217 


6 


HD 14489 


A2Ia f 


Allab 


Per OBI 


5.178±0.009 


0.369+0.004 


-0.110±0.000 


26/09/2001 


360 


246 


7 


HD20041 


AO la 


AO la 


Cam OBI 


5.795±0.019 


0.712±0.020 


0.090 


30/09/2001 


600 


234 


8 


HD21291 


B9 la f 


B9Ia 


Cam OBI 


4.213+0.019 


0.412±0.008 


-0.234+0.009 


26/09/2001 


2x240 


259 


9 


HD39970 


AO la 


AO la 


Field 


6.018+0.004 


0.386+0.005 


-0.192±0.010 


30/09/2001 


600 


270 


10 


HD46300 


AO lb' 


AO lb 


MonOBl 


4.498±0.008 


0.007+0.009 


-0.217±0.041 


29/09/2005 


180 


206 


11 


HD 186745 


B8Ia 


B8Ia 


VulOBl 


7.030±0.008 


0.930+0.002 


0.028+0.007 


25/09/2001 


900 


163 


12 


HD 187983 


Alia 


Alia 


Field 


5.590±0.026 


0.684±0.017 


0.173±0.149 


25/09/2001 


300 


187 


13 


HD 197345 


A2 Ia s 


A2 la 


Cyg OB7 


1.246±0.015 


0.092±0.007 


-0.233±0.008 


21/09/2005 


8x20 


798 


14 


HD202850 


B9 Iat/ 


B9Iab 


Cyg OB4 


4.233±0.009 


0.123+0.011 


-0.386±0.026 


29/09/2001 


120 


231 


15 


HD207260 


A2 lab 7 


A2Iab 


Cep OB 2 


4.289+0.007 


0.518+0.011 


0.119±0.018 


26/09/2001 


120 


370 


16 


HD207673 


A2 nV 


A2Ib 


Field 


6.467±0.005 


0.410±0.000 


0.060 


29/09/2001 


720 


195 


17 


HD208501 


B8HV 


B8Ib 


CepOB2 


5.796±0.004 


0.724±0.008 


-0.022+0.007 


26/09/2001 


480 


231 


18 


HD2 10221 


A3 lb 7 


A3 lb 


Field 


6.140±0.000 


0.414±0.017 


0.240±0.000 


26/09/2001 


720 


271 


19 


HD212593 


B9latf 


B9Ib 


Field 


4.569±0.018 


0.086+0.004 


-0.342±0.006 


29/09/2001 


2x180 


403 


20 


HD2 13470 


A3 la 


A3 lab 


Cep OBI 


6.650 


0.560 


0.240 


29/09/2001 


900 


249 


21 


BD+602582 


B8Iab 


B8Iab 


Cas OB2 


8.694+0.261 


0.770±0.016 


0.017±0.011 


27/09/2001 


2400 


140 


22 


HD223960 


A01a f 


B9Ia 


Cas OB5 


6.895±0.009 


0.715+0.009 


-0.050±0.047 


25/09/2001 


1200 


226 



Foces fl = 65 000 



23 HD195324 Al lb Allb Field 5.880+0.000 0.524+0.014 0.100 07/10/2001 2x1000 618 
Feros tf = 48 000 



24 


HD34085 


B8Ia« 


B8Ia 


OriOBl 


0.138+0.032 


-0.029±0.004 


-0.666+0.018 


14/11/1998 


20 


634 


25 


HD87737 


AO lb* 


AO lb 


Field 


3.486+0.053 


-0.026±0.015 


-0.206±0.028 


21/01/1999 


120 


440 


26 


HD91533 


A2 lab 


A2 lab 


Car OBI 


6.005±0.019 


0.318±0.011 


-0.075±0.034 


23/05/2007 


100 


229 


27 


HD111613 


A2 Iabe 


Alia 


CenOBl 


5.741±0.019 


0.384+0.022 


-0.088+0.026 


23/01/1999 


600 


376 


28 


HD 149076 


B8Iab 


B9Ib 


AraOBlb 


7.373±0.018 


0.485+0.009 


-0.118±0.021 


24/05/2007 


280 


230 


29 


HD 149077 


B9Ib 


AO lb 


AraOBla 


7.433±0.082 


0.470±0.021 


0.097±0.049 


24/05/2007 


310 


261 


30 


HD 165784 


A2/A3Iab 


A2Iab 


SgrOBl 


6.538±0.016 


0.856±0.005 


0.279±0.056 


09/07/2007 


140 


145 


31 


HD166167 


B9.5Ib 


AO lb 


SgrOBl 


8.605+0.009 


0.560+0.000 


0.036±0.047 


09/07/2007 


610 


117 



Uves tf = 80000 



32 


HD80057 


Allb 


Allab 


Vela OBI'' 


6.044±0.016 


0.285±0.006 


-0.117= 


=0.021 


24/02/2003 


2x50 


293 


33 


HD102878 


A2 lab 


A2Iab 


Cm OBI 


5.695±0.017 


0.265±0.009 


-0.119: 


=0.057 


06/01/2002 


54+139 


442 


34 


HD 105071 


B8 la/lab 


B8Iab 


Field 


6.316±0.024 


0.200±0.010 


-0.436= 


=0.034 


26/02/2002 


2x74 


382 


35 


HD 106068 


B8 la/lab 


B8Iab 


Field 


5.920±0.010 


0.297±0.007 


-0.284= 


=0.151 


20/07/2001 


2x53 


415 



Notes. <a) Blaha & Humphreys (1989); {b) Mermilliod & Mermilliod ( 1994); {c) note that the exposure times for UVES objects can vary in different 
wavelength bands; (< " adopte d from the SIMBAD data base at CDS, set in italics if a revisi on appears to b e re quired (see n ext column); (e) this 
work; (/) MK standards from lJohnson & Morgan! d!953h ; 00 anchor points of the MK system dGarrisonll 19941) ; (/l) lReed j2000h . 



reduction procedure a median filter was applied to the raw data 
to remove cosmic ray hits without losing spectral information. 
Special attention had to be paid to the merging of orders, be- 
cause key features like the hydrogen lines extend over adjacent 
orders and their analysis is sensitive to line shape. Moreover, 
meticulous normalisation procedures were applied (within the 
reduction pipeline), since imprecise continuum determination 
may lead to incorrect abundance and parameter determination. 
Comparison of our data with longslit spectra of fundamental 
MK s tandards obtained within the NStars project dGrav et al.l 
l200l was possible for some cases and showed good agree- 
ment. Figure[TJexemplifies such a comparison of a NStars project 
1.8 A-resolution spectrum around the B aimer lines Hy to H8 of 
Deneb (taken at the Dark Sky Observatory) with our Foces spec- 



2 http : //stellar . phys . appstate . edu/ 



trum, which was artificially degraded in resolution to match that 
of the NStars project spectrum. Finally, for several objects mul- 
tiple exposures were taken and combined to enhance the quality 
and filter out artifacts. 

Additional objects were observed in 2007 at the European 
Southern Observatory on La Silla, u sing Feros (Fiber -fed 
Extended Range Optical Spectrograph, iKaufer et al.l 1 19991) on 
the 2.2 m telescope. The spectra cover the wavelength range 
from -3600 to 9200 A at 48 000. The raw spectra were 
mainly processed by the ESO automatic reduction pipeline. In 
addition to spectra of our science targets, a template spectrum of 
a subdwarf B star was obtained and reduced with the pipeline. 
By dividing this spectrum through a well-fitting model we could 
successfully filter out all artifacts from the data reduction pro- 
cess, providing - after proper smoothing to remove localised 
noise relics - a template function for the continuum rectifica- 
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4320 4330 4340 4350 4360 




0.2 b , , , : 

4080 4090 4100 4110 4120 




0-2 b , , , : 

3950 3960 3970 3980 3990 




0-2 1 



3870 3880 3890 3900 3910 

A (A) 

Fig. 1. Comparison of our Foces spectrum (full line) with a 
longslit spectrum from the NStars project dGrav et alj [2003, 
dots) in the region of the Balmer lines Hy to H8 for Deneb 
(A2 la) to assess the quality of the normalisation of our echelle 
spectra. Note that the Foces spectrum was was artificially de- 
graded in resolution to match the 1.8 A resolution of the longslit 
data. 



tion of the science target spectra. This procedure provides an 
objective and reproducible means for the normalisation process. 

Additional high-quality spectra of bright BA-type super- 
giants were extracte d from the database of the UVES Paranal 
Observatory Project dBagnulo et al.ll2003l) . These data combine 
a high resolution of R w 80 000 and high S/N with large wave- 
length coverage from 3040 to 10400 A, completing our sample. 

All together, this comprises the most comprehensive sample 
of Galactic BA-type supergiants^ with high-quality spectra con- 
sidered for quantitative analysis to date. 

The high-quality spectra allow a much closer look to be 
taken on trends in spectral line strengths and line ratios than 
possible with traditional classification spectra at much lower res- 
olution. We used the opportunity to reassess the spectral clas- 
sification of the sample stars. Starting point of our approach 
were the anchor points of the MK system (as identified in 
Table [Q, supplemented b y MK primary standards as given by 
iJohnson & Morgan! (Il953l) . see also Table [T] These cover about 
half of our sample stars. Spectral types for the remaining stars 
were assessed on basis of the helium and metal lines, luminos- 
ity classes by the width of the Balmer lines. In general, good 
agreement with the classification as obtained from SIMBAD was 
found. Maximum changes indicated by our inspection, for about 
half of the non-MK stars, amount to one spectral subtype or to 
one subtype (higher or lower) within the supergiant luminosity 
class. Our reclassification is indicated in Table [T] We note that 

3 Besides a supergiant nature also a scenario as a p ost-AGB candidate 
is considered for HD 195324 (Szczerba et al. 2007) b ecause of its in- 
frared excess measured by IRAS (laschek et al. 1991; Oudmaiie r et al] 
119921) . We will address the evolutionary state of this interesting object 
based on surface abundance patterns in Paper II. 



two of the primary MK standards (HD 14489 and HD 223960) 
were found to differ significantly in spectral morphology from 
the stars of similar original spectral type. Also the luminosity ap- 
pears overestimated for HD 212593, which shows a nearly sym- 
metric Ha absorption profile. We therefore propose a (slight) 
reclassification for these stars as well, based on the available 
high-resolution spectra. 

Spectra of several sample stars are displayed in Fig. [2] for 
three strategic wavelength regions around Hy, Ha and Pal3. A 
sequence in spectral type from B8 to A3 is shown, the earliest 
and latest types observed here. Note the rapid strengthening of 
the metal lines towards lower effective temperatures, and their 
enormous increase in number. The spectra were shifted to the 
laboratory rest frame and the major spectral features are identi- 
fied, many of which were analysed in our project. 

Not unexpected are the asymmetric Ha line profiles, indicat- 
ing the presence of a stellar wind of varying strength in the dif- 
ferent stars. P-Cygni-like profiles like in /JOri are an exception, 
occurring only in the highest-luminosity objects (being most ex- 
treme in HD 12953). Typical for many stars of the sample is ex- 
cess emission in the red wing of Ha like for the AO- A3 stars in 
Fig. [2] Only the luminosity class lb stars show Ha completely 
in absorpti on, with only small as ymmetries being present. See 
the atlas of Verdug o et al.l ( fl9~9 9a) for more examples on A-type 
supergiant spectra. 

In addition to the spectra, which constitute the princi- 
pal data for the analysis, various (spectro-)photometric data 
were adopted from the literature for constructing spectral 
energ y distributions. Johnson t/BV-magnitudes were taken 
from [Mermilliod & Mermilliodl d 19941) (see Table Q), which 
are means of previously published photoelectric data, and 
7//A"-magnitudes from the Two M icron All Sky Survey 
dCutri et alj|2003l: ISkrutskieetal] [2006, 2MASS). Additionally, 
flux-calibrated, low-dispersion spectra observed with the 
International Ultraviolet Explorer (IUE) were extracted from 
the MAST archiveQ, where available (i.e. for 16 objects). For 
seven objects only high-resolution IUE spectra are present in 
the MAST archive, which are nevertheless useful for our pur- 
poses as they were observed using a large aperture (i.e. flux- 
losses should have not occured). These data cover the range from 
1 150 to 1980 A for the short (SW) and from 1850 to 3290 A for 
the long wavelength (LW) range camera. Typically, both wave- 
length ranges were observed the same day. A summary of the 
individual spectra used in the present work (data ID number 
and observation date) is given in Table [2] Finally, photometric 
data in the Stromgren syste m was adopted from the catalog of 
lHauck & M ermilliodl d!998l not tabulated here). 

3. Quantitative analysis 

3.1. Models and analysis methodology 

Advanced models and robust analysis techniques are required 
in order to achieve the targeted accuracy and precision in stel- 
lar parameter and abundance det ermination. A suitable analy- 
sis methodology was presented by Przybilla et al.l (I2006I) . which 
we adopted for the present work, with some extensions as sum- 
marised below. The method is based on a hybrid non-LTE ap- 
proach, in which non-LTE line-formation computations are per- 
formed on top of LTE model atmospheres^ 

4 http://archive.stsci.edu/ 

5 First results obtained with hydrodynami c, fully line-blanketed non- 
LTE atmospheres computed with Cmfgen ( Chesneau et al. 2010) and 
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8690 
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Fig. 2. Normalised spectra of an T e ff -sequence of sample stars around Hy (left), Ha (middle) and Pa 13 (right panel), from the hot 
(top) to the cool end of the parameter range investigated here (bottom). The major spectral features are identified, short vertical 
marks indicate Tin lines. The Ho- region is contaminated by narrow telluric lines. Vertical shifts of the spectra by multiples of 0.5 
have been applied for clarity. 



Table 2. IUE spectra used in this study. 



Table 3. Model atoms for non-LTE calculations. 



# 


Object 


SW 


Date 


LW 


Date 


1 


HD12301 


P07282 


01/12/1979 


R06276 


01/12/1979 


2 


HD12953 


P42698 


12/10/1991 


P21249 


15/09/1991 


3 


HD 13476 






P29471" 


07/11/1994 


4 


HD 13744 






P31590" 


14/10/1995 


5 


HD 14433 






P09166" 


24/09/1986 


6 


HD 14489 


P21812 


19/12/1983 






7 


HD20041 


P56064 


09/10/1995 


P31581 


09/10/1995 


8 


HD21291 


P07280 


01/12/1979 


R06274 


01/12/1979 


9 


HD39970 


P56171 


09/11/1995 


P31675 


09/11/1995 


10 


HD46300 


P56165 


08/11/1995 


P31667 


08/11/1995 


12 


HD 187983 


P48688 


19/09/1993 


P26414 


19/09/1993 


13 


HD197345 


P09133 


26/05/1980 


R07864 


26/05/1980 


14 


HD202850 


PI 5099 


25/09/1981 


R11614 


23/09/1981 


15 


HD207260 


P03368 


17/11/1978 


R02957 


17/11/1978 


16 


HD207673 






P31682" 


10/11/1995 


17 


HD208501 


P55805 


03/09/1995 


P31403 


03/09/1995 


18 


HD2 10221 


P18682 


28/11/1982 


R14745 


28/11/1982 


19 


HD2 12593 


P33852 


03/07/1988 


P13556 


03/07/1988 


20 


HD2 13470 






P29702" 


13/12/1994 


22 


HD223960 






R08994" 


11/10/1980 


24 


HD34085 


P31880 


18/09/1987 


PI 1654 


18/09/1987 


25 


HD87737 


P08566 


26/03/1980 


R07305 


26/03/1980 


33 


HD 102878 






P28097" 


09/05/1994 



Notes. ia) high-resolution spectrum 

In brief, LTE model atm ospheres were computed using the 
Atlas9 code (lKuruczll993ah with some additional modifications 
necessary to allow model convergence close to the Eddington 



Ion 



Model atom 



H 

He i 

Ci/n 

Ni/n 

Oi/n 

Mg ifu 

S n/m 

Tin 

Fen 



Przvbilla & Butler ( 
Przybilla (2005) 

Przvbilla et al. (2001b) : Nieva & Przv billa (2006. 2008) 
Przvbilla & Butler (2001) 

Przvbilla et al. (2000); Becker & Butler (1988) . updated 
Przvbilla et al. (200 1 a) 
Vranckcn ct al. (1996). updated 
Becker (1998) 
Becker (1998) 



PoWR (W.-R. Hamann, priv. comm.) indicate that the photospheric 
spectra of BA-type supergiants are described well by our approach. 



limit (Przvbilla et al. 2001b). Line blanketing w as accounted fo r 
via Kurucz' Opacity Distribution Functions dKuruczl fl9 93b). 
Based on these model atmospheres the non-LTE calculations 
were performed with recent versions of Detail and Surface 
(Giddines 11981b iButler & Giddingslll985l both updated by K. 
Butler). The coupled radiative transfer and statistical equilibrium 
equations were solved with Detail, providing non-LTE level 
populations. Subsequently, synthetic spectra were computed 
with Surface, using refined line-broadening theories. State-of- 
the-art model atoms according to Table [3] were adopted for the 
calculation of the model grids. 

Overall, around 25 000 model spectra per element were com- 
bined into 5-dimensional grids. The parameter space covered 
ranges from 8300 K to 15 500 K in effective temperature T e g 
(in steps of 250-500 K), from 2.50 in surface gravity logg (cgs 
units) to the convergence limit at 0.95 (lower r e ff-limit) to 1.90 
(upper r e ff-limit, in steps of 0. 1 dex), from 3 km s to 8 km s 
in microturbulence £ (lkms -1 steps) and from 0.09 to 0.15 in 
helium abundance (by number, steps of 0.015). Finally, suffi- 
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Fig. 3. Example for a set of diagnostic lines used in the determination of the atmospheric parameters T e g an d log g in the supergiant 
HD 106068 (B8 lab), consisting of hydrogen, nitrogen and oxygen lines. The models calculated for our final parameters (red lines) 
are compared to observation (black lines). The same values of oxygen and nitrogen abundance are adopted for all lines, regardless 
of ionization stage. The ionization equilibria and the Balmer lines unaffected by the stellar wind are reproduced simultaneously. 
Note that various line blends of stellar, interstellar and terrestrial nature seen in the plots are exempted from the fitting process. 



cient coverage in elemental abundance (typically 1 dex in steps 
of 0.25 dex) for the respective lines investigated is provided. 

Models of Doppler-shift profiles were computed by disk in- 
tegration, using the methods described in iGravl d2005l) . We di- 
vided the apparent disk of a star in 16 million parts of equal 



size. For each disk position, the radial-tangential velocity dis- 
tribution, representing macroturbulenc^ was projected on the 



A physical explanation of the 'macroturbulent' velocity parameter 
was only recently suggested for hot stars, lik ely being the collective ef - 
fect of pulsations dAerts et ail2009h . see also Simon- Diaz et al.l (2010). 
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Fig. 4. Abundance values derived from single lines are pre- 
sented versus respective equivalent widths for a selection of lines 
from neutral and ionized nitrogen and iron in the spectrum of 
HD 46300 (AO lb). The lcr-scatter around the mean value is in- 
dicated by the gray bands. 

line of sight, assuming an equal amount of radial and tangen- 
tial motion. The result was then Doppler-shifted according to 
stellar rotation. Finally, we applied a linear limb-darkening law 
to determine the weighting factors for integration. The resulting 
profiles are characterised by three parameters: the projected ro- 
tational velocity v sin i, the macroturbulent velocity and the 
limb-darkening coefficient e. 

Theoretical and observed line profiles were compared via 
the software p ackage Spas (Spectral Plotting and Analysis Suite, 
iHirschl 120091) . It provides the means to interpolate between 
model grid points for up to three parameters simultaneously and 
allows to apply instrumental and macrobroadening functions to 
the resulting theoretical profiles. Furthermore, the program uses 
the downhill simplex algorithm to minimise x 1 in order to find a 
good fit to the observed spectrum. 

3.2. Spectroscopic indicators 

In order to find a globally satisfying solution, it was necessary 
to derive the basic atmospheric parameters T s g, log g, £ and y in 
an iterative process, utilising different spectral indicators on the 
way. Here we discuss the individual steps of this procedure. 

3.2.1 . Effective temperature and surface gravity 

The hydrogen lines and the ionization equilibria of the metals 
react sensitive to variations of effective temperature and sur- 
face gravity. This dependency was demonstrated for the Stark- 
broaden ed hydrogen Balmer l ines as well as for e.g. the Mg i/ ii 
lines by IPrzvbilla eTail (I2006I) and lSchiller & Przvbillal d2008). 
While neither the hydrogen lines nor a single ionization equilib- 
rium alone are sufficient to constrain T e ff and log g unambigu- 
ously due to a degeneracy in the solutions, a combination of at 
least two of them is, for constant £ and y. Note that in the low- 
temperature regime of our sample, where N n lines are too weak 
to be analysed, the C i/n ionization equilibrium is often available 
to confirm the results obtained from the Mg i/ii lines, while in 
the high-temperature regime O i/n becomes available in addition 
to Ni/ii. For the majority of the stars two or in some cases even 
three metal ionization equilibria can thus be utilised for the at- 
mospheric parameter determination. He (see Fig. [2j and, in the 
most luminous objects, even H/3 and Hy can be affected by the 
presence of a significant stellar wind. Consequently, these lines 
have to be omitted in such cases for the derivation of basic at- 
mospheric parameters in our hydrostatic approach. H/3 and Hy 



are therefore considered for most stars, and several of the higher 
Balmer lines are analysed in all stars, with the exact number de- 
pending on the wavelength coverage of the spectrograph, and the 
S/N reached in the blue orders. If ionization equilibria for more 
than one element are available for analysis, the results show typ- 
ically high consistency. An example for the agreement of a final 
solution with observation is shown in Fig. [3] 

3.2.2. Micro-turbulence and helium abundance 

The microturbulent velocity is usually determined by demand- 
ing that the abundances indicated by the lines within an ion are 
independent of equivalent width. This is basically equivalent to 
finding the microturbulence value that minimises the line-to-line 
scatter in abundance, which is our preferred method (as we use 
line-profile fits for this like for the rest of the analysis). Within 
the framework of the present paper the nitrogen spectra are the 
most useful indicators for microturbulence as they provide many 
lines of widely differing strength, while e.g. the more extensive 
Fe ii line analysis (to be discussed in Paper III), among others, 
can be used to verify the ^-determination. An example is shown 
in Fig. [4] using the traditional form of illustration. Equivalent 
widths were determined using direct integration over the line 
profile. Linear regression curves to the data points are consistent 
with a slope of zero in both cases. A single value for £ was found 
to be consistent for all ele ments analysed in non-L TE in the indi- 
vidual stars, see Fig. 13 of IPrzvbilla etail (120061) for a compre- 
hensive discussion in the case of HD 87737. Note that while the 
B aimer-lines are virtually insensitive to microturbulence varia- 
tions, this quantity influences the parameter determination based 
on ionization equilibria. 

A change in the atmospheric helium abundance gives rise to 
a modified mean molecular weight of the atmospheric plasma, 
which has simi lar effects on model predictions than a variation 
of logg dKudritzki|[l973t IPrzvbilla et all 12006). Therefore, the 
helium abundance y has to be constrained simultaneously with 
the other atmospheric parameters. This is done via line-profile 
fitting of helium lines. Results for the individual Hei line abun- 
dances in the sample stars can be found in Paper II. 

In practice, estimates for £ and y are used to determine initial 
values for T e ff and log g in a first iteration step (using elemental 
abundances as third fit parameter), which are in turn adopted to 
derive f and y in a next step. Convergence is quickly achieved 
using our comprehensive model grids and Spas. 

3.2.3. Projected rotational velocity and macroturbulence 

In order to perform line-profile-fitting it is necessary that the 
models are able to reproduce the line-shapes accurately. To 
achieve this our model spectra are convolved with the functions 
for rotational and macroturbulent broadening computed as de- 
scribed earlier. A value of 0.5 for the limb-darkening coefficient 
e was adopted throughout. Tests showed little sensitivity for this 
parameter within acceptable bounds, which is only slowly vary- 
ing around our canonical value throughout the atmospheric pa- 
rameter and spectral wav elength range investigated here (see e.g. 
IWade&Rucinskilll985l) . 

We chose the same set of highly suitable metal lines to de- 
rive vsini and ( in most of our targets, including the Mgn 
line at 4390 A, the Mgi line at 5183 A, the Fen pair around 
6148 A, and the O i triplet around 6157 A. However, for the two 
hottest objects in our sample only the Sn line at 5354 A was 
used. Examples from the fitting are shown in the upper panels 
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Fig.5. Examples for the derivation of vsin/ and £ from the Oi triplet around 6157 A (left) in HD 207673 and the Fen pair around 
6148 A (right) in HD 14433. The upper panels show contour plots of a goodness-of-fit parameter in (vsin/, £) space, lower val- 
ues representing a better fit. The lower panels display the best line-profile-fits with and without macroturbulence (left and right, 
respectively). 



of Fig. [5] displaying a goodness-of-fit parameter in (v sin i, £) 
space. While in theory this should allow us to identify the opti- 
mum values for v sin i and f , we found disentanglement of the ef- 
fects from rotational velocity and macroturbulence to be difficult 
because of the rather smeared contours. This is similar to other 
recent studies of rotational broadening in BA-type supergiants 
(IVerdugo et al.lll999bl:lRvans et al]l2002l) . Values of v sin i and f 
representing the best fits show a scatter from line to line because 
of this degeneracy in the solutions. Fits without macroturbulence 
contribution agree very well in the resulting v sin i for different 
absorption features. However, they can not reproduce observa- 
tions in such a satisfactory way like fits that consider both ef- 
fects, see the lower panels of Fig. [5] for a comparison. Overall, 
the uncertainties of vsin / and f are typically about +5kms _1 
because of the degeneracy of the solutions, and about ±3 km s _1 
for the stars with the sharpest lines. 

3.2.4. Abundances and metallicity 

Abundances for all ionization stages involved in the atmo- 
spheric parameter determination were derived, as summarised 
in TablelU For this, the individual lines of an ion were fitted and 
the mean of the resulting abundances was taken. For the sake of 
brevity we postpone a presentation of individual line abundances 
for the sample stars to Papers II and EH, where the complete data 
will be discussed, irrespective of the availability of ionization 
equilibria. 

Overall, good to excellent agreement between different ion- 
ization stages of the elements was found, even in cases where 
more than one ionization equilibrium was evaluated, indicating 
that major bias is likely absent in our atmospheric parameter 
determination. The statistical uncertainties arising from noise 
prove to be small due to the high quality of our spectra. Only 
for the weakest lines used in our analysis they are comparable 
to the line-to-line scatter. Therefore we chose the lcr standard 



deviation based on the individual line abundances as a conserva- 
tive error estimate for the ion abundances. In this approach it is 
difficult to assign an uncertainty value for an ion where only one 
line is available for analysis. An inspection of the data in Table|4] 
indicates that a characteristic estimate in such a case may be 
O.lOdex. 

Low star-to-star scatter in metallicity was found in re- 
cent analyses of early B-type dwarfs - the progenitors of BA- 
type super giants on the main se q uence - in the solar neigh- 
bourhood (IPrzvbilla et all 120081: iNieva & Simon-Dfazl 1201 lb 
iNieva & Przvbillal 1201 2|) . using similar analysis techniques. 
Given the l ow sensitivity of ou r analysis to small changes in 
metallicity (Przvbilla et al. 2006), this parameter was held fixed 
throughout our model grid to slightly subsolar values (with re- 
spect to the solar standard of lGrevesse & Sauval 1998) in accor- 
dance with the B-dwarf results (the cosmic abundance standard). 
This approximation was proved correct a-posteriori by our abun- 
dance analysis for most of the sample stars. We reanalysed the 
most metal-rich stars of our sample later on by means of fine- 
tuned micro-grids, adopting a higher value for the metallicity in 
order to rule out any possible bias due to this. 

3.2.5. Spectral energy distribution 

As a verification of our spectroscopic analysis we compared the 
Atlas9 fluxes computed for our final parameters with photomet- 
ric data in the optical and near-IR, and UV spectrophotometry. 
Thus, we investigated whether the models reproduce also the 
spectral energy distributions (SEDs) of the stars, and as a conse- 
quence also their global energy output. 

We transformed the UBV and JHK magnitudes into absolute 
fluxes using zeropoin ts from [ Bessell et all (Il998l) for Johnson 
photometry and from lCohen et all (120031) for the 2MASS pho- 
tometric system. Note that we lowered the U photometric zero 
point by 3% based on discrepancies between theory and observa- 



8 



M. Fimstein and N. Przybilla: Quantitative Spectroscopy of Galactic BA-type Supergiants. I 




3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 3.0 3.2 3.4 3.5 3.8 4.0 4.2 4.4 

log X (A) log X (A) 



Fig. 6. Examples of comparisons of Atlas9 model fluxes (black lines) with UV-spectrophotometry from the IUE satellite (gray lines) 
and with photometric measurements - UBV from Mermilliod & Mermilliod (1994), JHK from 2MASS - for two of our sample 
stars. The SEDs shown are dereddened according to the values of E(B - V) and Ry in Table|U and normalised in V. 



tion that were found bv lBohlin & Gillilandl (2004) in the respec- 
tive part of the UV-flux for the photometric standard star Vega. 

The observed fluxes needed to be dereddened for the com- 
parison with the model flux es. For this, we adop ted the mean 
interstellar extinction law o f lCardelhetaTl (ll989). which is de- 
pending on one parameter, the ratio of total-to-selective extinc- 
tion Ry = Ay I E(B - V). The colour excess E(B - V) was deter- 
mined per star as the difference between the theoretical colour 
computed from the Atlas9 flux and the observed colour. Then, 
Ry was varied until a good overall match of the theoretical and 
the observed SED was achieved, see Fig.|6]for examples of final 
results. The typical uncertainty in Ry is +0.2. This could be done 
also for the twelve stars without IUE spectrophotometry, though 
lacking the very thorough constraints by the UV fluxes. 

In general, as good agreement between the model and the ob- 
served SED is found as exemplified in Fig. [6] Notable deviations 
occur only in three cases. An IR excess is found for HD 21291, 
while the UV-flux is in excellent agreement. We attribute this 
to HD21291 being part of a c lose visual binary separated by 
only 2.39" dPrieur et al.l 120081) . which could not be resolved 
by 2MASS. An UV excess in HD 187983 can be explained 
by adding flux from an early B-type main-sequence compan- 
io n, which is con sistent with a radial-velocity variation found 
bv lHendrvl ( ll981l) . There remain some problems in reproducing 
the SED of HD 12953, the most luminous object in our sample, 
simila r to those found for HD 92207 (AOIae) by iPrzvbilla et al l 
(2006). Further investigations are required to address the issue, 
but these are beyond the scope of the present work. 

Finally, a by-product of our investigations are bolometric 
corrections B.C. for the sample stars. They were derived from 
the fluxes of our final models, see Sect. I5.4l for a discussion. 

4. Results 

4.1. Atmospheric parameters of the sample stars 

Our results on the atmospheric parameters and some derived 
properties of the sample stars are summarised in Table |4] 
For each object our internal sequential number and its HD- 
designation are given, followed by our determinations of ef- 
fective temperature, surface gravity, microturbulent, projected 
rotational and macroturbulent velocity, colour excess, ratio of 
total-to-selective extinction and bolometric correction. Finally, 
surface abundances of helium and - for the cases where the re- 
spective ionization equilibria were employed for the atmospheric 
parameter determination - abundances of metal ions are sum- 



marised. The error margins denote lcr-uncertainties (an estimate 
of O.lOdex may be applied for ion abundances where only one 
line could be analysed, i.e. those data without uncertainty indica- 
tor). Uncertainies in vsin/ and £ amount typically to +5kms~ 1 
(+3kms -1 for the stars with the sharpest lines) because of the 
degeneracy of the solution in the fitting of vsin; and see 
Sect. 13.2.31 and +0.2 in Ry, see Sect. 13.2.51 The uncertainties of 
the bolometric corrections vary from ~0.02mag for the coolest 
sample stars to ~0.05mag for the hottest stars, see Sect. 15. 41 for 
a more detailed discussion. This is the most comprehensive col- 
lection of such data from a homogeneous non-LTE analysis of 
Galactic BA-type supergiants to date. As expected, a continuous 
distribution over the temperature range from 8400 to 12 700 K is 
obtained, displaying no signs of systematic shifts or gaps stem- 
ming from the use of different temperature indicators. 

Tests with a bootstrapping method implemented in Spas in- 
dicate formal uncertainties of less than 0.05 dex in logg and 1% 
in T e g, when all other parameters are held constant at the val- 
ues of our final solution. This applies even for the lowest-quality 
spectra of our sample (at S/N -ratio of about 120). Fits to individ- 
ual Balmer lines - the main indicators for \ogg throughout the 
sample - show a scatter of less than 0.05 dex around the mean 
value in this parameter in the individual stars. This suggests that 
major systematic effects from order merging and normalisation 
in the hydrogen lines may be absent. 

The line-to-line scatter in abundances derived from individ- 
ual metal lines, however, indicates low-level systematic uncer- 
tainties from within our model atoms and the normalisation of 
weak lines. Small selection effects may also emerge in cases 
where only a few lines in one ion of the two analysed for es- 
tablishing ionization balance are available. Therefore, when line 
sets of sufficient quality for two or more elements were available, 
we compared different solutions using only one of the available 
ionization equilibria and the Balmer lines. While still showing 
good agreement, the results imply slightly higher uncertainties 
than discussed earlier. Hence we chose our uncertainty estimates 
carefully across the parameter space, remaining conservative in 
the values given in Table|4] The largest discrepancies were found 
for the case when the last residual Mg i lines (Mg i A5 1 83 A be- 
ing the strongest) are about to vanish. In consequence, the largest 
uncertainties occur in the transition region when N i/n starts to 
replace the Mg i/n-balance. Overall, the typical conservative un- 
certainties in r e ft and logg are therefore ~2% (1 50-200 K) and 
0. 10 dex (i.e. ~25% in gravity), respectively, accounting for both 
statistical and systematic errors. 
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Table 4. Stellar parameters of the program stars. 
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Notes. (o) Typical uncertainties are ±5 kms and ±3 kms 1 for the stars with the sharpest lines. ib) Typical uncertainties are ±0.2. <c) Typical uncertainties vary from ~0.02mag for the coolest 
sample stars to ~0.05 mag for the hottest stars. {d) In cases where only one spectral line could be analysed for an ion (no error margin given), a good estimate for the uncertainty is O.lOdex. 
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Table 5. Comparison of our atmospheric parameters (r e ff 
logg) for objects in common with previous studies. 
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Object 


Takeda" 


Venn'' 


McErlean r 


this work 
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HD12301 


12500, 2.30 




14000, 2.15 


12600, 2.15 
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9000, 1.50 


8400, 1.20 




8500, 1.40 
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HD 14489 




9000, 1.40 




9350, 1.45 
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12000, 1.80 




11500, 1.60 


10800, 1.65 
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HD39970 


10500, 1.70 






10300, 1.70 


10 


HD46300 


10000, 2.00 


9700, 2.10 




10000, 2.15 


13 
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9500, 1.50 






8700, 1.20 
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11500, 1.80 






10800, 1.85 
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9500, 1.50 






8800, 1.35 


16 
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9300, 1.75 




9250, 1.80 


17 
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12500, 2.30 




13000, 1.80 


12700, 1.85 
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8200, 1.30 




8400, 1.40 
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11200, 2.10 


23 
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9200, 1.85 


24 
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13000, 1.75 


12100, 1.75 


25 
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10200, 1.90 


9700, 2.00 
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Spectra 
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Fig. 7. Comparison of our results for the individual sample stars 
with reference spectral type-r e ff scales. Luminosity classes and 
reference scales are encoded according to the legend. A typical 
error bar is indicated to the lower left. See Sect. I5.1l for a discus- 
sion. 



Our accuracy in microturbulence determination is limited by 
the stepsize of 1 kms -1 in our model grid. The values derived 
from different elemental species are consistent within this er- 
ror. The derived helium abundances vary between 10 and 14% 
in number fraction, with errors of around 20% in these values. 
Note, however, that the high sensitivity of the helium lines to 
r e ff-variations makes the values prone to substantial systematic 
shifts even from small errors in the parameter determination. 
Finally, metal abundances are determined on average with about 
10-20% uncertainty. Abundances from lines of the major ioniza- 
tion stages are usually insensitive to variations of and log g 
within the given uncertainties. 

4.2. Comparison with previous analyses 

Samples of about 10-20 Galactic supergiants in the BA-star 
regim e were subj ect to non-LTE analyses in three previous stud- 
ies (IVennlll995a[ IMcErlean et al.lll999HTakeda & Takada-Hidail 
2000). We list their effective temperatures and surface gravities 
for stars in common with the present work in Table [5] Typical 
uncertainties in T eft /logg are 200K/0.2dex, 500-1000 K/0.5 dex 
and 1000 K/0.2 d ex in these studies, respectively. 

IVennl <Q995a) used a similar method for the parameter de- 
termination in A-type supergiants, i.e., utilising the Balmer- and 
Mg i/n-lines, however without a simultaneous derivation of the 
He abundance in her hybrid non-LTE approach. While there are 
some differences, good overall agreement is found in and 
log g within the uncertaintie s. 

Takeda & Takada-Hidai (2000) tried to construct a T e g/ log g 
vs. spectral type relation, utilising published data. They empha- 
sised, however, that this approach is subject to considerable un- 
certainties. Indeed, we find notable deviations between our val- 
ues for both T e ff and logg and their data, though agreement 
is still found when considering their rather large uncertainty 
ranges. 

Finally , there is some overlap with the B-type supergiants 
studied by McErlean et al. ( 1999). Their estimates for logg, in- 
ferred from the Balmer lines, are in excellent agreement with 
ours. Yet, we find a systematic shift in T e ff towards higher val- 
ues compared to our results. Their temperatures are based on the 
non-LTE Si n/m ionization balance for the stars in common (ex- 



cept for HD21291, where photometry was used as substitute). 
However, the non-LTE model atmospheres used for their analy- 
sis did not consider the effects of metal line-blanketing, which 
we identify as the likely cause of the r e ff-differences. 



5. First applications 

5. 1 . Spectral type-T en relation 

Empirical spectral-type-r e (j relations provide important start- 
ing points for all kinds of stellar studies and for quantitative 
spectroscopy in particular, and they are therefore an essential 
part of the reference literature on stellar properties. Our high- 
precision/high-accuracy dataset facilitates to reassess the exist- 
ing knowledge in the BA-type supergiant regime in view of im- 
proved models and analysis techniques. 

Effective temperatures of our sample stars (see Table 2]) as a 
function of their spectral type (according to the refined classifi- 
cation presented in Table [TJ are displa yed in Fig.[7l In compari- 
son w ith established reference work (ISchmidt-Kalerl[l982l : ICoxl 
2000) we find a significantly steeper spectral-type-r e ff relation, 
i.e. higher T e ff at spectral types B8 and B9, and lower T e s at A2 
and A3. As no apparent correlation of r e ff with luminosity class 
(LC) is indicated in Fig. [7] we compute average r e jf-values from 
all stars of a spectral subtype to provide refined reference values 
which are presented in Table|6] 

With typically five to seven objects per spectral type it would 
be desirable to increase the sample size in order to get a more ro- 
bust determination, in particular of the uncertainties. The distri- 
bution of the data per spectral bin is not Gaussian at the moment. 
Note that the formal lcr-uncertainties for spectral types Al and 
A3 in Table|6]are lower than the uncertainties from the individual 
star analyses. There is little overlap of r e ff-values for individual 
stars of different spectral types. Overall it appears that our re- 
fined spectral-type- T e g relation allows effective temperatures to 
be estimated with a typical lcr-uncertainty of about 2-3%. 

Obviously, the work has to be extended in the future to pro- 
vide a unified pict ure accounting for recent developments con- 
cerning hotter (e.g. Rep olust et alJl2Q04l iMarkova & Pulsl[2 008) 
and cooler Galactic supergiants (e.g. lLevesque et aUl2005h . 
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Table 6. Spectral type-r e ff-(B - V)o— (2> - y)o scales & B.C. values for Galactic BA-type supergiants. 



reference values from literature this work'' 



Sp.T. 


1 eff 


(B - V) a 


1 eff 


(B - V) " 


0>-yV 


B.C. a ' Ba 




(B 


-V) 


(b-y) 


B.C. 




K 


mag 


K 


mag 


mag 


mag 


K 




mag 


mag 


mag 


B6 


13000 


-0.08 








-0.76 












B7 


12200 


-0.05 








-0.66 












B8 


11200 


-0.03 


11100 


-0.03 




-0.54 


12200±410 


-0.09: 


tO.Ol 


0.001±0.009 


-0.66±0.08 


B9 


10300 


-0.02 








-0.40 


10920+220 


-0.06: 


t0.02 


0.011+0.012 


-0.31 ±0.04 


AO 


9730 


-0.01 


9980 


-0.01 


0.032±0.019 


-0.29 


9840±290 


-0.05: 


tO.Ol 


0.017±0.013 


-0.19±0.05 


Al 


9230 


0.02 








-0.20 


9240± 80 


-0.01: 


t0.02 


0.043±0.014 


-0.09±0.03 


A2 


9080 


0.03 


9380 


0.03 


0.051 ±0.011 


-0.16 


8960+200 


0.00: 


t0.02 


0.051±0.015 


-0.03±0.05 


A3 


8770 


0.06 






0.074±0.021 


-0.09 


8430+ 60 


0.01: 


tO.Ol 


0.056±0.005 


+0.08±0.02 


A5 


8510 


0.09 


8610 


0.09 


0.062+0.017 


-0.01 













Notes. (D) Schmidt-Kaler 1 19821) : m ICoxl d2OO0h : (c) lGravl il992h . for LC lb; (rf) corrected by +0 ra 12 to match the zeropoint of the present work; 
(c) see Sects. [57TT45.4l for a discussion. 
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Schmidt-Kaler (1982) 
Cox (2000) 
Flower (1996) 
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HD12953 



X 


la 


A 


lab 


• 


lb 



Groy (1992) 



(B-V) 

Fig. 8. Comparison of observational results for the individual 
sample stars with established reference (B - V)o-r e ff relations 
from the literature, according to the legend. A typical error bar 
is indicated. See Sect. l5.2l for a discussion. 



(b-y)o 

Fig. 9. Comparison of observational results for the sample stars 
with a reference (b-y)o-T e g relation from the literature, accord 
ing to the legend. A typical error bar is indicated. See Sect 
for a discussion. 



5.2. Ce-Vjo-Tgff scale 

Closely related to the spectral-type-r e ff relation is the question 
of the behaviour of intrinsic Johnson colour (B - V)o with spec- 
tral type, or more precisely, with effective temperature. We de- 
termined (B - V)o from the observed (B-V) colour (see TableQ]), 
corrected for the colour excess E(B - V) according to Table [4] 

Our effective temperatures as a function of (B - V)o are dis- 
play ed in Fig. HQ together with establi shed reference relations 
fromlSchmidt-Kalerl d 19821 for LC Iabl. lFlowerl d 19961 LC I) and 



ICoxl d200d LCD~Our data are systematically bluer by (F03- 
0^06 than indicated by the reference relationships. The super- 
giants become bluer with decreasing luminosity (i.e. smaller ra- 
dius) for a given r e ff, the decrement being about O. n 01 per lumi- 
nosity subclass. As we do not have enough objects per spectral 
type to compute meaningful means for the individual luminosity 
subclasses, we provide (B - V)o-values averaged over the entire 
LC I in Tabled Note that the outlier HD 12953 - the by far the 
most luminous star of our sample - has been excluded from this. 

The loci of the three luminosity subclasses of BA-type 
supergiants in Fig. [8] may be approximated by a quadratic fit 
function, yielding the following relations 

(B-V)o(mag) = 

la: 50.668-24.585 log r eff + 2.977 (log r e ff) 2 (1) 



lab : 34.368 - 16.529 logr eff + 1.981 (logr eff ) 2 
lb : 49.565-24.114 log T eff + 2.927 (log T eS ) 2 



(2) 
(3) 



with an area of validity of about 3.92 < log T e g < 4.10. Like for 
the case of the spectral type-r e ff relation it would be desirable 
to extend this to hotter and cooler supergiants. 

5.3. (b-y)o-T e ff scale 



Our effective temperatures as a function of Stromgren (b-y )n are 
displa yed in Fig. [9] together with a reference relation from Grav 
(1992, for LC lb). There is an excellent match of Gray's relation 
with our data. Note the drastic change of the slope of the (b—y)o- 
r e ff-relation for bluer (b-y)o. In general, the supergiants become 
bluer with decreasing luminosity for a given r e ff, the decrement 
being about OT'Ol per luminosity subclass. Like for the previous 
case of the Johnson data, we provide (b - y)o-values averaged 
over the entire LCI in Table [6] The outlier HD 12953, which 
appears exceptionally blue in (b - y)o, has been excluded. 

The loci of the three luminosity subclasses of BA-type 
supergiants in Fig. [9] may be approximated by a quadratic fit 
function, yielding the following relations 

(b - y) (mag) = 
la : 30.222 - 14.633 log r eff + 1 .77 1 (log r eff ) 2 



(4) 
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Fig. 10. Comparison of bolometric corrections determined here 
with reference relations from the literature, according to the leg- 
end. The error bar denotes a typical uncertainty in r e ff and a 
maximum uncertainty in B.C.. See Sect. l5.4l for a discussion. 

lab: 12.147- 5.669 log r eff + 0.659 (log T eS ) 2 (5) 
lb: 42.265-20.736 log r eff + 2.543 (log r eff ) 2 (6) 

with an area of validity of about 3.92 < log r e ff < 4.10. 



The derived values are compared to the reference calibra- 
tions of Schmidt- Kaler ( 1982, for stars of luminosity class lab), 
which are reproduced by ICoxl ([2000), and of iFlowen ((1996, 
for supergiants, his Table 4) in Fig. [TO] In order to adjust for 
the different zeropoints used, we correct the data of Schmidt- 
Kaler/Cox byO^lzfl. The resulting differences can reach values 
up to 0™13. A general trend towards larger B.C-values is seen in 
our data, except for the hottest objects in the sample, which have 
slightly lower bolometric corrections. The match of the Flower 
relation (adjusted by +0^01 because of the different zeropoints) 
on the other hand is good, with small deviations becoming ap- 
parent at the highest and lowest temperatures considered here. 

As no significant deviations occur due to luminosity clas- 
sification, except possibly for our most luminous sample star 
HD 12953, we average over all members of a spectral subtype to 
determine improved B.C. reference values. Our values are com- 
pared to the reference data from the literature in Table [6] 

5.5. Photometric atmospheric parameter determination 

As one often wishes to analyse larger samples of stars, easy-to- 
apply and fast stellar parameter indicators are in high demand. 
Photometric r e fr -indicators are among the preferred ones, be- 
cause of the easy accessibility of photometric data. We discuss 
empirical relations for the two most common photometric sys- 
tems, for Johnson and for Stromgren photometry. The latter also 
allows surface gravities to be estimated photometrically. 



5.4. Bolometric corrections 

We determined bolometric corrections from the model fluxes 
for ea ch individual star. Following the approach of lBessell et al.l 
(1998), the model flux at the V-band effective wavelength was 
converted into a Johnson V magnitude in the usual way by adopt- 
ing a zeropoint flux = 3.636 x 10~ 20 erg cmT 2 s _1 Hz (for 
Vega). The integral over the total model spectral energy distribu- 
tion yielded a bolometric magnitude (assuming an absolute so- 
lar bolometric magnitude M^ ol = +4.74). The difference between 
the bolometric and the V magnitudes calculated in this way then 
provided the bolometric correction, which can be expressed as 



B.C. = M bol -M V = C- 101og(r eff /r° ff ) + 2.5 logff v , 



(7) 



where C ■■ 

r° ff = 5777 K, and H. 



12.854 for a solar radius fl = 6.95 x 10 lu cm and 
is the Eddington flux at the effective V- 
band frequency. This way, the correct empirical B.C. for the 
standard star Vega is reproduced, which lies exactly in the r e ff- 
range spanned by the sample stars. Note that for this reason we 
have not renormalised the bolometric corrections to the mini- 
mum value of zero in order to avoid positiv e Zj.C.-values for 
some late-A and F-type stars, as suggested by iBuser & Kurucz] 



1978), se e also the dis cussion of this topic by iBessell et alJ 



1998b and lTorresi d2010h . As a result, our B.C. values are about 



0™12 larger than the ones provided earlier by us when intro- 
ducin g the modelling techniques and analysis methodology used 
here (Przybill a et al.l l2006h . The B.C-values determined from 
Eq. |7J summarised in Table |4] deviate less than 0^03 from 
the recent analytical fit formula of Kudritz ki et al.l (12008b. their 
Eq. 6), which have used the same approach. The uncertainties 
of the bolometric corrections are dominated by the r e ff-errors. 
For the coolest sample stars, i.e. close to the minimum of the 
r e ff-5.C.-curve, the uncertainties for individual objects amount 
to ~0.02mag, while they reach up to ~0.05mag for maximum 
r e ff-uncertainties in the hottest stars. 



The reddening-free Johnson Q-index as T eff -indicator. 

Effective temperature c alibrations based on the reddening-free 
g-index (|Johnsonll958l) . Q = (U-B)- X(B - V) with X = E(U- 
B)/E(B - V), have come into wide use recently, as Johnson pho- 
tometry is available for most stars. Examples in our context en- 
compass studies of the evolutionary progenitors of BA-type su- 
pergiants, OB-type stars on th e main sequence (e.g. lDaflon et al.l 
1999; Lvubim kov et al.| | 2002[). and their coole r siblings, late-A 
to G-type supergiants (ILvubimko v et al.ll2010|). Usually, a stan- 
dard value of 0.72 is adopted for X, but see Johnson! (Q/958). 

We test the applicability of the Q-index as possible 
r e ff-indicator for BA-type supergiants in Fig. QT| Our spec- 
troscopically derived r e jf-values are displayed as a function 
of Q. A trend with luminosity subclass becomes apparent, the 
more luminous la objects being coolest and the lb objects being 
hottest at a given Q. The loci of the three luminosity subclasses 
in the diagram may be approximated by a quadratic fit function, 
yielding the following relations 



7f ff (10 3 K) 



la 
lab 
lb 



8.691 +3.447g+ 13.843 Q 2 , 
8.435 + 3.430 2+ 17.672 Q 2 , 
8.725 -0.026 2+ 12.631 Q 2 , 



- 0.30 > Q > -0.65 (8) 
-0.15 > Q > -0.60 (9) 

- 0.05 > Q > -0.55 (10) 



with their respective area of validity. The outlier HD 12953 was 
excluded from the fitting procedure. For clarity, only the fit func- 
tion for lab supergiants is visualised in Fig. QT] At first glance, 
rather tight relations are found. 

The differences between the spectroscopically derived ef- 
fective temperatures T^ c and the g-based computed with 



7 ISchmidt-Kaieil dl982T) ad opts M°, = -0™19, while the value pro- 
vided by IBessell et al.l dl998h is -0™07. Note that the Schmidt-Kale r 
data is subject to a further internal inconsistency by 0" 1 02 (Torres 2010), 
which would improve the comparison with our values only slightly. 
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Fig. 11. Relation between the reddening-free Johnson Q-index, 
Q = (U - B) - 0.7 '2(B - V), and the spectroscopic r e ff-values of 
the sample stars. A typical error bar is indicated. The dashed line 
represents the regression line of our £)-based r e ff-calibration for 
lab supergiants. 



o 
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Fig. 12. Upper panel: comparison of our spectroscopically de- 
rived r^ ec with the Q-based computed with Eqs.lBffTOl The 
dashed line indicates the 1 : 1 relation. The error bar shows a typ- 
ical conservative uncertainty from our spectroscopic determina- 
tion, and the uncertainty in T® s , accounting for typical uncertain- 
ties in the colours only. Lower panel: percent difference of the 
two r e ff-values for the individual objects. These quantify the sys- 
tematic differences from application of our Q-T e g calibrations. 
The dotted lines indicate the lcr-scatter range. 



Eqs. [BTTTOl are quantified in Fig. Q~2] The sample stars follow 
the 1:1 relation rather tightly, with the lcr-scatter of the com- 
puted around the spectroscopic reference values amounting 
to more than 3%. This is the systematic uncertainty adherent to 
an application of the method. In addition, random errors because 
of the uncertainties in the colours also have to be considered, 
typically amounting to more than twice the lcr-uncertainty of 
the spectroscopically derived values, i.e. another ~5% (see the 
error bar in the upper panel of Fig. [12] and note the different 
scale projection of the axes). In consequence, the resulting error 
margins render this simple approach of Q-based effective tem- 
peratures not competitive with spectroscopic determinations for 



X 


a 




A 


ab 




• 


b 





a\ 



HD12953 X 



4< v A 
X 



0.0 0.2 0.4 0.6 



1.0 1.2 



Fig. 13. As Fig. [TTJ but for the relation between the reddening- 
free Stromgren [c;]-index, [ci] - c\ - 0.20(b - y) where c\ - 
(u - v) - (v - b)), and the spectroscopic r e ff-values of the sample 
stars. 
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Fig. 14. As Fig. Q~2] but for the comparison of our spectroscop- 
ically derived r^ ec with the [ci]-based T^' computed with 



Eqs. [TTJ4J1 



BA-type supergiants. Moreover, better starting points for an it- 
erative refinement of the stellar parameters are already obtained 
from the spectral type-r e ff relation established in Sect. 15.11 

The reddening-free Stromgren [ci]-index as T e tf-indicator. 

Data for 31 o f our sample superg ia nts w ere found in the 
catalogue of lHauck & Mermilliodl (1 19981) on Stromgren 
uvbyB photometry, facilitating the reddening-free [cj]-index 
[ci] = c\ -0.20(£> —y) to be computed, with c\ — (u - v) - (v-b). 
Figure [13] visualises the sample stars in the [ci]-r e ff plane. The 
loci of the three luminosity subclasses may be approximated 
by a quadratic fit function as before, yielding the following 
relations 



la 
lab 
lb 



r' e ff l] (io 3 K) = 

14.610- 11.706 [ci] +5.855 [ci] 2 ,0.15 < [a] < 0.90(11) 
15.741 - 12.710 [ci] +5.536 [ci] 2 ,0.35 < [a] < 1.10(12) 
14.310- 5.893 [ci] + 1.092 [a] 2 ,0.25 < [a] < 1.25(13) 
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Fig. 15. Relation between Stromgren B and the spectroscopically 
derived flux-weighted gravities loggp of the sample stars. The 
error bar in B is an estimate based on the restricted information 
available from the Stromgren data. The dashed line represents 
the linear regression to the data, ignoring the outlier HD 210221. 

with their respective area of validity. For clarity, only the fit func- 
tion for lab supergiants is visualised in Fig.Q~3] The clear outlier 
HD 12953 was ignored in the determination of the fit coefficients 
in Eq.QT] The Stromgren [ci]-r e ff relations are obviously tighter 
than the Johnson Q-T^ relations. In fact, a practical common 
relation for BA-type supergiants of luminosity class la and lab 
can be given with only a negligible loss in accuracy for the range 
0.15 < [ci] < 1.10: 

(10 3 K) = 14.613 - 10.565 [ci] + 4.523 [ci] 2 . (14) 

In analogy to Fig. [12] the differences between r^ ec and 

the [cj]-based 7™ computed with Eqs. [jHT - [T3"1 are quantified in 

Fig. [14] The lcr-scatter of the computed Tzj! around the spec- 
troscopic reference values amounts to about 2% only in this case 
(systematic uncertainty of the method). The random errors due 
to Stromgren colour uncertainties are also drastically reduced 
compared to Johnson photometry, amounting to about 1%. We 
conclude that the reddening-free [c\\-index provides highly use- 
ful starting values for T e g in high-precision/accuracy analyses 
of all except the most luminous BA-type supergiants. 

As the reddening-free Stromgren [ci]-index measures 
the Balmer decrement at the temperatures of the sam- 
ple stars, it could act as a sub stitute for the Balmer 
jump method (Kudritzki et al. 2008b) based on flux-calibrated 
low/intermediate spectroscopy, that has been used for studies 
of BA-type supergiants in galaxies beyond the Local Group re- 
cently. Extended exposures in a spectral region with usually 
low instrumental sensitivity (decrease of mirror coating reflec- 
tivity and CCD sensitivity towards the optical-UV) could thus 
be avoided. However, practical challenges to be solved first are 
the determination of the metallicity dependency of the [c\]-T e g 
relation and the availability of Stromgren filter sets at large tele- 
scopes. 

Stromgren B as a photometric gravity indicator. The 

Stromgren /Hndex is a poor direct indicator for surface grav- 

8 Note that the Balmer jump shows some sensitivity to variations of 
surface gravity as well. 



Fig. 16. Upper panel: comparison of our spectroscopically de- 
rived log gp Pec with the /3-based log g^ computed with Eq. [16] 
The dashed line indicates the 1 : 1 relation. The error bar shows 
a typical uncertainty from our spectroscopic determination, and 
the uncertainty in log g^., accounting for typical uncertainties in 
B only. Lower panel: difference of the two log gp-values for the 
individual objects, in dex. These quantify the systematic differ- 
ences from application of our B-log gp calibration. The dotted 
lines indicate the lcr-scatter range. 

ity in the parameter range of BA-type supergiants because the 
strength of the Balmer lines depends not only on \ogg - pres- 
sure broadening via the linear Stark effect - alone but also on 
r e (j - via (de-)population of the n — 2 level of neutral hydrogen 
by excitation/ionization. On the other hand, there is a good cor- 
relation with the flux-weighted gravity 

log g F = log g - 4 log T eSA , (15) 

as illustrated by Fig. [BJ Here, T eSA = r eff /10000K. The 
flux -weighted gravity is a mea sure for the stellar luminos- 
ity ([Kudritzki et alJl2003l l2008bl) . with lower values indicating 
higher luminosity. Linear regression of the data yields 

log^ = -12.4548 + 5.449 B. (16) 

One statistical outlier, HD 210221, was excluded from the fit 
based on its distance to the remainder of the data points. The 
random lcr-uncertainty due to errors in the measurements of 
B is typically 0.03 dex, and the application of Eq. [16] yields 
the flux-weighted gravity to within 0.1 3 dex (systematic lcr- 
uncertainties), see Fig. [TBI 

The conversion of log g^ to the /3-based surface gravity log g 13 
is facilitated by correction with the [ci] -based effective temper- 
ature, as established before. The uncertainties of T]%r contribute 
little to the error budget of log g@, resulting in error margins of 
+0.04 dex (lcrrandom)+0.13dex(lcr systematic uncertainty). 

In particular the [ci]-r e ff (assisted by a determination of the lu- 
minosity subclass from spectroscopy) and the y6-log gp relations 
derived here should facilitate to provide useful starting values of 
r e ff and log g for detailed quantitative analyses of normal BA- 
type supergiants. Exceptions may be the most luminous objects 
(those with pronounced Ha P-Cygni profiles like HD 12953), 
which as a class of their own fall out of relations given here. 

The photometric starting values must be refined - preferen- 
tially through an iterative methodology like the one discussed 
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Fig. 17. The sample stars with Stromgren photometric data in the 
T e s—logg plane. Solutions from the spectroscopic determination 
are shown as dots, open circles mark values that are obtained 
using photometric indicators according to our calibrations using 
[ci] and ft. The two corresponding solutions are interconnected. 
Typical error bars are indicated in the upper left. Dotted lines 
mark evolution track s for rotating stars at metallicity Z = 0.014 
dEkstrom et alJl2012l) . from 9 M to 32 M (bottom to top). The 
9 M model shows a blue loop, partially displayed here. 



in Sect. [3]- in the spectroscopic analysis if high accuracy and 
precision is desired in the derivation of all dependent quantities. 
The reason for this is illustrated in Fig. [17] which shows both the 
results of the atmospheric parameter determination from our de- 
tailed spectroscopic analysis and from application of the photo- 
metric indicators discussed here in the r e ff-log g plane. Despite 
"inconspicuous" error bars (see Fig.fTTt. the photometric indica- 
tors without a subsequent refinement would yield different solu- 
tions for the fundamental stellar parameters (mass, radius, lumi- 
nosity) and for elemental abundances. Symptoms of the presence 
of such systematics would become evident in any more detailed 
analysis, as a failure to establish the ionization balance for multi- 
ple elements simultaneously, and very likely a mismatch of the- 
oretical and observed profiles for some hydrogen lines. Results 
biased in such a way can easily lead to misinterpretations and 
incorrect conclusio ns in the astrophysical context, see Fig. 5 of 
Sim on-Dfazl (1201 Oh for a good example on the effect for the el- 
emental abundance determination, and his discussion of this. In 
terms of fundamental stellar parameters, we see e.g. that evolu- 
tionary masses from the comparison with stellar evolution tracks 
in Fig.[T7]could be biased by up to ~30%. 



5.6. Tracers for studies of the ISM 

Properties of the interstellar medium (ISM) are usually inves- 
tigated by using early-type stars as background light sources, 
facilitating to study absorption and scattering by the interven- 
ing interstellar material along the line of sight. While the topic 
is too broad for the scope of the present paper, we nevertheless 
want to point out the usefulness of BA-type supergiants for such 
studies briefly, in particular for the example of the ratio of total- 
to-selective extinction Ry. 

Our derived 7?y-values for the sample stars (see TableUJi clus- 
ter around the frequentl y quoted standard value of Ry — 3.1 (e.g. 
ISavage & Mathis 1979), with an rms scatter of 0.3. The data are 
displayed in Fig. [18] as a function of position in the Galactic 
plane (Galactic coordinates I and b are used, note that two 
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Fig. 18. The distribution of our sample supergiants in the 
Galactic plane, coded according to the legend for values of the 
total-to-selective extinction ratio Ry along the line-of-sight. See 
Sect.l5.6lfor a discussion. 



stars fall outside the displayed range, HD 34085 and HD 87737, 
which are located at higher Galactic latitude). We find the objects 
with higher values than Ry = 3.1 concentrated in the southern 
Milky Way around 260° to 360° Galactic longitude, and many 
objects with lower values in the opposite direction, at /»60 
to 110° . This is in good agreement with previous findings of 
IWhittell(ll977h . 

While our sample is too small to provide any significant ex - 
tension to modern studies (see e.g. Win klerl 19971: IWegnerl2 003). 
which consider hundreds of OB stars, we want to draw the at- 
tention to the potential of BA-type supergiants for such inves- 
tigations. Despite their rarity, they are highly useful for prob- 
ing lines-of-sight towards very distant (or highly reddened) ob- 
jects in the Milky Way because of their much higher intrinsic 
visual brightness. Moreover, their usefulness for simple investi- 
gations like present ed here extends to distances even far beyon d 
the Milky Way (e.g. lKudritzki et alJ2008bll20TllU et al.l2 009). 
It is reassuring in this context to find the ^-values for many 
independent Galactic sightlines to cluster close to the canonical 
value. A ratio of total-to-selective extinction of 3.1 may there- 
fore be a good assumption for many extragalactic environments, 
where full SED information of individual stars, required for the 
determination of Ry, is unavailable at present. 

In addition, luminous BA-type supergiants facilitate also 
more sophisticated studies of the ISM in other galaxies 
when high/intermediate-resolution spectroscopy becomes feasi- 
ble. Examples are investigations of the neutral interstellar gas 
via the Na D lines (for B-type supergiants, with no stellar contri- 
bution t o these lines) or of di ffuse interstellar absorption bands 
(see e.g. lCordiner et alJ[2008allbh . 



6. Summary 

A sample of 35 bright Galactic BA-type supergiants was in- 
troduced. Quantitative non-LTE analyses of high-resolution and 
high-S/N spectra were performed with the aim to provide a 
homogeneous set of atmospheric parameters at highest accu- 
racy and precision. The study provides the most comprehensive 
dataset on effective temperatures, surface gravities, helium abun- 
dances, microturbulent, macroturbulent and rotational velocities 
of Galactic BA-type supergiants so far. In addition, the inter- 
stellar reddening and the ratio of total-to-selective extinction to- 
wards the sample stars were determined. 
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First applications of the data show that established rela- 
tions from the reference literature for BA-type supergiants are 
outdated. An improved empirical spectral-type-r e ff and T e g- 
B.C. scale was derived, as well as new calibrations of intrinsic 
Johnson (B - V)o and Stromgren (b —y)o colours as a function of 
effective temperature were provided. It would be desirable to ex- 
tend such studies to supergiants of earlier as well as later spectral 
types, with overall improved number statistics. 

Photometric r e ff-determinations based on the reddening-free 
Johnson Q-index were found to be of limited use for high- 
precision/accuracy studies of BA-type supergiants because of 
large errors of typically ±5%(lcr statistical)+3% (lcr system- 
atic), compared to a spectroscopically achieved precision of 1- 
2% (combined statistical and systematic uncertainty with our 
methodology). On the other hand, the reddening-free Stromgren 
[ci]-index and Stromgren B are highly promising for deriving 
good starting points for further quantitative investigations, with 
uncertainties of + 1% +2.5% in T e tf, and ±0.04±0.13dex in \ogg 
(lcr-statistical, lcr-systematic, respectively). Finally, the poten- 
tial of BA-type supergiants as tools for studying ISM properties, 
also in other galaxies, was briefly addressed. 

Besides the immediate objectives, this paper provides the 
preparatory work for addressing important questions of mod- 
ern astrophysics. Tight observational constraints on stellar evo- 
lution in form of fundamental stellar parameters and light ele- 
ment abundances as tracers of rotational mixing will be derived 
in Paper II, and constraints on Galactic chemical evolution, via 
abundance gradients of the thin disk, in Paper III. 
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